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I. Introduction 



This report describes the development of coursework in automatic controls 
for the Mechanical Engineering (ME) Department which took place during FY 85 
as a continuation of work begun in FY 84 [Ref* 1]. This work was motivated by 
the Educational Skill Requirement (ESR) levied by NAVSEA as follows: 

"[Teach] basic understanding of automatic control systems and 

their application to Integrated Propulsion Plant Control." 

The ME Dept, approach to the ESR was to identify two required ME courses 
to teach the basics. In addition, the department began to investigate 
elective courses which could support an automatic controls thesis specialty 
option. Further, the total offering was to be consistent with a widespread 
practice of controls teaching. This report is organized into six sections: 

The present introductory section is followed by a brief summary of previous 
work; the third and fourth sections discuss FY 85 course and laboratory 
development, respectively; the fifth section presents related research; and 
the last section contains conclusions. 

II. Previous Work 

Work accomplished in FY 84 led to the following conclusions : 

1. Two required courses in controls are needed to satisfy the ESR. 

There are so many basic concepts in systems and controls which need 
to be introduced in these courses that some understanding of 
applications will only be achieved through a period of study 
equivalent to approximately two quarter-courses in length. 

2. A specialization level of understanding (4XXX, graduate level 
course) is not appropriate to either of the required courses, as per 
the ESR. 
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3. ME faculty should teach both required courses. Automatic control is 
a fundamental area of mechanical engineering and, as such, all of the 
ME faculty in the area of design should be able to teach the two 
courses • 

4. Some "hands-on" controller investigation by the students is desirable 
to complete the students’ understanding of basic concepts. 

5. A design-oriented approach to the second required course should be 

taken in order to keep the students in touch with the goal of control 
design: a successfully controlled system (e.g. a controlled power 

plant in accordance with specifications). 

Based on the preceding conclusions, the following courses were 
recommended : 

Required 

ME3801 Linear Automatic Control. 

(to replace EC3413, Fudamentals of Automatic Control) 

Introduction to linear controller analysis methods. Classical 
methods for single-input-single-output (S1S0) plants. 

ME3802 Marine Control Systems. 

(to replace ME4802, Marine Propulsion Control) 

Introduction to controller design practice. Linearization of 
nonlinear plants. Linear controllers for multiple-input 
multiple-output (MIMO) plants. Study of Integrated Throttle 
Control on the DD963 and other developmental MIMO controllers. 

Elective 

ME4801 Fluid Power Control. 

Analysis and design of hydraulically actuated control sy terns. 

ME4803 Advanced Topics in Controls. 

System Identification (parameter identification). Nonlinear 
controller design. Design of microprocessor-based controllers. 

ME4902 Reading Course in Controls. Selected, specialized topics for 
individual study with the professor. 
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In addition, it was recommended that thesis topics in controls be 
identified within the areas of marine propulsion, servomechanisms, and 
robotics . 

III. FY85 Course Development 

ME3801, Linear Automatic Control, is being offered for the first time 
during AY 86-1. The course syllabus is included as Appendix I. The objective 
of the course is to provide the students with the basic analytical tools 
necessary to evaluate a controlled system composed of a linear controller and 
a linear, SISO plant. The course is taught from a fundamental viewpoint with 
general application studies. 

ME4802, Marine Propulsion Control, is the course which is now taught to 
satisfy the ESR regarding integrated control. This course has evolved into a 
design-oriented course, compared to ME3801 which is analysis-oriented, with 

4 

special emphasis on control of marine propulsion plants. The course was 
taught twice during AY 85 and the latest course outline is included in this 
report as Appendix II. 

The study of the design problem in ME4802 was begun by examining Navy 
Controls Specifications, both the general specifications for ships of the U.S. 
Navy and the Proposed DDG51 design specif ications . A DD963 marine gas turbine 
propulsion plant was then evaluated to identify cause and effect 
relationships, plant inputs and outputs, and control variables. The 
implementation of Integrated Throttle Control on the DD963 was next discussed 
in order to illustrate a contemporary control strategy implementation (a 
recent research report was used to do this, "The Naval Gas Turbine Ship 
Propulsion Dynamics and Control Systems Research and Development Program," 
SNAME Trans., Vol. 90, 1982, pp. 321-338^. Following this specific 
introductory material, a more general approach was taken to the remainder of 
the course. 
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The course approach for ME4802 was developed to address the modeling and 
control design tasks in the order that they must be faced by a controls 
designer. In this way, it was felt that the students can learn what must be 
done, when it is done, and how it is done by studying the controller design 
process. Specifically, the following design tasks were discussed: 

1. Specification for control design. 

2. Evaluation of plant function. 

3. Plant mathematical modeling. 

4. Plant model validation - open loop simulation and experimentation. 

5. Selection of control strategy. 

6. Selection of actuators and sensors. 

7. Dynamic modeling of actuators and sensors. 

8. Selection of controller action. 

9. Theoretical controller design. 

10. Controller validation - closed loop simulation and comparison to 
specifications . 

11. Prototype. 

The design process was repeated four times for four different propulsion 
plants. The course emphasized tasks 1, 2, 5, 8, 9 and 10 for the various 
controller/piant combinations. Plant models were supplied in the interest of 
expediency (tasks 3 and 4), as were actuator and sensor models (tasks 6 
and 7). The propulsion plant was assumed to be given and unchangeable, thus 
allowing the class to fully concentrate on the controller design process. 

The theory content of the course was developed around the idea that the 
plant is usually nonlinear, but that it can be linearized. This is a good 
approach for marine propulsion systems. Consequently, linear controller 
design methods were taught based on plant input/output classification. For 
singie-input-single-output (SISO) plants, PID and compensation methods were 
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discussed. For multiple-input-multiple-output (MIMO) plants, transfer matrix 
compensators and linear quadiatic regulators were discussed. 

The course was concluded with a discussion of the implementation of 
digital controllers. The following hardware topics were discussed: analog- 

to-digital conversion of sampled data; sampling rate effects; 
computer/controller recursion formulas; and digital-to-analog conversion of 
controller signals • 

The present "laboratories" for ME4802 are really ten homework projects, 
each of which consists of general theory homework problems from the course 
text plus one or more controller design tasks for a marine gas turbine 
propulsion plant. During the quarter, each student completed four different 
controller designs for the plant. In order to give the students a feeling of 
reality, the design tasks were related to the NPS marine gas turbine 
propulsion emulator (ref. 2). Copies of the projects are included as Appendix 
III. Further work needs to be done to support these homework projects with 
true lab projects and "hands on" hardware involvement (more about this 
below) . 

The elective course ME4801, Fluid Power Control, has been successfully 
offered in the past the only change has been in lab development which will be 
discussed below. 

ME 4512, Advanced Dynamics, has not been taught for some time at the NPS. 
However, a growing student interest in robotics at the school has caused a 
renewed student interest in the course material. The course is now under 
review and will be offered in the winter quarter, AY86-2. 
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ME4803, Advanced Topics in Controls, is in the formative stage. It could 



be offered as student interest indicates and as lab equipment becomes 
available • 

IV. Laboratory Development 

Lab development is proceeding in conjunction with the NAVSEA program for 
lab hardware improvement now being conducted at the NPS. There are two phases 
to this program: 



Phase I. 


Hardware was purchased to support the fluid power control course. 


ME4801. 


The hardware included: 


1 . 


One Hewlett Packard HP85 controller and data acquisition system. 


2. 


A fluid power bench for synthesizing hydraulic control experiment 
setups . 


3. 


Measurement instrumentation for measuring fluid power variables. 


4. 


An armdraulic table-top instruction robot. 



Supporting this hardware purchase, a series of experiments has been 
designed to illustrate some important features of hydraulic controls 
operation. Five experiments and their objectives were defined as follows 
(for more details, see Appendix IV) . 



i. 


Hydraulic actuator static performance. To determine the static 
performance of a typical rotary actuator (pump driven motor). 


2. 


Flow control valve static performance. To characterize a flow 
control valve of typical configuration (e.g. four-way, three-end, 
critical center or open center spool valve). 


3. 


Dynamic performance of a hydraulic power element, valve controlled 
position. To evaluate the performance of a hydraulic power element, 
to determine its performance characteristics (leakage coefficients, 
hydraulic spring rate, damping ratio, etc,) and to compare these 
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results with theory, 

4. Dynamic performance of a position-control servomechanism. To 
examine the dynamic performance of the servo and to evaluate the 
effects of various design options. 

5. Dynamic performance of a velocity-control servomechanism. To 
examine the dynamic performance of the servo and to evaluate the 
effects of various design options. 

Phase II . This phase of the NAVSEA/NPS Lab hardware improvement program is 
now in the proposal stage. Hardware items for experiments to support ME4802 
(soon to be ME3802, Marine Control Systems) and ME4803, Advanced Controls 
Topics, have been proposed. The following experiments are envisioned: 

A. For ME4802, Marine Propulsion Control: 

1. Diesel control lab. This will be a table-top experiment using a 
microprocessor-based diesel engine simulator and a programmable 
controller to investigate transfer matrix compensation. 

2. Steam engine control lab. Another table-top apparatus will be used 
to simulate steam engine operation and to investigate the use of a 
linear quadratic regulator for control. 

3. Gas turbine control lab. This experiment will investigate the 
implementation issues of sampling rate, A/D conversion, D/A 
conversion, and computer algorithms. A table-top 

microprocessor-based simulation of the marine gas turbine needs to 
be developed to conduct this experimentation. 

B. For ME4803, Advanced Topics in Control, the distinguishing feature between 
these labs and those for ME4802 will be the emphasis on controller design. 
The labs for ME4802 will be more instructor prepared, "cookbook'* fashion, 
and will allow the students to concentrate on the input-output 
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relationships between the controller and the plant for various controller 
types. The labs for ME4803 will require the students to investigate the 
implementation issues in controller design. Two types of labs will be 
developed : 

1. Microprocessor-based controller lab. To develop an understanding ‘of 
the power and limitations computer/controllers. 

2. Analog controller lab. To study the power and limitations of analog 
controllers • 

V. Related Research 

ME Department research in marine propulsion control has been looking into 
issues of control of diesel and marine gas turbine systems. The following 
theses have been advised during FY85: 

1. "Marine Propulsion Load Simulation," P. N. Johnson, MSME Thesis, June 
1985. 

2. "Modeling of Marine Gas Turbine Components," J. Roger, MSME Thesis, 
December 1985 (expected). 

3. "Modern Control of a Marine Gas Turbine," V. Herda , ME Engineer’s 
Thesis, June 1986 (expected). 

4. "System Identification and Control of a Internal Combustion Engine," 
T. Violette, MSME Thesis, December 1985 (expected). 

ME Department research in robotics has been very active in response to 
the Naval Surface Weapons Center (NSWC/White Oak) Robotics Lab. A long-term 
program with NSWC is now being developed to ensure a continuing source of 
thesis topics for NPS students. Much of the work so far has been stimulated 
by the NSWC firefighter robot project. Robotics research projects during FY85 
included the following: 
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1. "Optimal Control of Robotic Mechanisms," D. L. Smith, NPS Foundation 
Research • 

2. "Firefighter Robot Prototype Development," D. L. Smith, NSWC proposed 
research for FY 86. 

3. "Linearized Controller Design for a Revoiute Robot," D. Lewis, MSME 
Thesis, December 1985 (expected). 

4. "Load Measurements for a Firefighter Tool," R. Yobs, MSME Thesis, 
March 1986 (expected). 

5. "Simulation of High Speed Dynamics for a Rigid Revolute Robot," W. 
McCarthy, MSME Thesis, December 1985 (expected). 

6. "Modeling of Flexible Link Dynamics," R. Petroka, ME Engineer 1 s 
Thesis, June 1986 (expected). 

VI. Conclusions 

1. The required ME coursework in controls is becoming a well-integrated 
two course sequence which will produce students who have a good basic 
understanding of automatic control systems and their application to 
propulsion plant control. 

2. Coherent elective coursework is developing in such a way that it 
offers to students the opportunity to prepare for a worthwhile thesis 
research project based on an ME controls coursework specialization. 

3. Promising research topics in marine propulsion control and robotics 
have begun to open up. Student interest in these problem areas seems 
to be quite strong, as does sponsoring lab interest. 
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ME 3801 - LINEAR AUTOMATIC CONTROL 
TENTATIVE SCHEDULE — FALL, 1985 



Week 

of 



Subject Matter 



Assi gned 
Read 1 ng * 



Assi gned 
Probl ems 
(Cat. B) 



1 


30 


Sep . 


Introduction. Mathemat i - 

cal background. Laplace 

Transf orms. 


Preface, 
Ch. 1, 2. 


-L 1 


o 


7 


Oct. 


Math, background continued. 
Math, models of physical 
systems . 


Ch. 2, 4. 

(Ch. 3) 


2—4 — 2—9 . 
4-1 - 4.4. 


3 


14 


Oct. 


Math, models continued. 


Ch. 4. 
(4-6) 


4-5 - 4-7, 
4-9 - 4.12, 
4-14. 


4 


21 


Oct. 


Quiz. Controllers and 

basic control actions. 


Ch . 5 • 


5-1 - 5-5. 


5 


28 


Oct. 


Basic control actions 
conti nued . 


Ch. 5. 
(5-6) 


5-6 - 5-9. 


6 


4 


Nov . 


Transient response 
anal ys.i s. 


Ch. 6. 
(6-7) 


6-1, 6-3 - 
6-9. 


7 


1 1 


Nov. 


Quiz. Transient response 
continued. Error 
anal ysi s. 


Ch . 6 , 7 . 

-.7-4) 


6- 10 - 6-13, 

7- 1 - 7-4. 


a 


18 


Nov. 


Root locus methods. 


Ch . a . 


8-1 - 9-3, 
8-5, 8-8, 3-1 


9 


26 


Nov . 


Frequency response methods. 


Ch. 9. 
(9-4) 


9-1 - 9-3, 9- 
9-6. 


10 




Dec . 


Quiz. Frequency response 

methods conti nued . 


Ch . 9 . 


9-8 - 9-11, 
9-13. 


1 1 


q 


Dec . 


Design and compensat 1 on . 


Ch. 10. 





Text: Ogata, "Modern Control Engineering," Prentice-Hal 1 , 1970 

Chapters and sections in parentheses are recommended. 

Hours: Tu, Th , Fri 1310-1400, Wed 1510-1700 (all in Ha 109). 
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APPENDIX II 




NAVAL POSTGRADUATE SCHOOL 
Monterey, California 



ME 4802: 


Marine Propulsion Control (3-2) 


Prerequisite : 


EE 3413, Fundamentals of Automatic Control 


Instructor : 


Professor David Smith 
HA 207A, X3383 


Office Hours: 


To be announced 


Course 

Pescri ption: 


Propulsion systems overview. 

Modeling and simulation of propulsion performance. 
Control systems design and implementation. 

Case studies of the DD 963 and DDG51. 


Course 
Objectives : 


To familiarize students with the control design process. 




To introduce students to modern control methods for 
multi variable Vsy stems . 




To demonstrate the use of classical and modern control 
strategies for marine propulsion plants. 


Text : 


Modern Control Engineering, by K. Ogata, Prentice Hall Inc., 
1970. Course notes 


Homework : 


Assigned to aid understanding of concepts. 

Will be collected. 

Solutions will be posted on the second deck of Halligan Hall. 


Labs : 


Design - oriented lab projects have been developed to allow 
students to practice the full control design process. Written 
lab reports will be submitted weekly according to the course 
outline. Lab reports are due at the beginning of the lab 
period. Selected students will orally present either homework 
or lab work during the scheduled lab period. 


Course 

Grading: 


Pass /Fai 1 

Course grade will be based upon lab performance. 
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Course Outline 



ME 4802 
Summer , 
Week 

1 



2 



3 



4 



5 



Marine Propulsion Conrol 



1985 

Date 

7/8 



7/15 



7/22 



7/29 



8/5 



8/12 



Topic 



• Introduction 

• Control Specifications 

• Plant Function: 

Multiport Analysis 

Lab - Lab report format 

• SISO Plant Modeling: 

Linearization 
T.F. from data 
Simulation, CSMP 



Ref 

notes 

Ch.l, 10.1 
notes 



11.1, 4.1, 4.3 
9.9, notes 
notes 



Lab - Project 1: Multiport Analysis and Linearization 



• Classical Control Design: 
Control Action 
PID Tuning 
Compensation 



5.1, 5.5 
notes 

10 . 1 , 10 . 2 , 10.6 



Lab - Project 1 : Open Loop Plant Model Validation 

• Modern Control Design 14.1 

• MIM0 Plant Modelling: - 

State Space Modelling 14.2 

State Space Modelling notes 



Lab - Project 1 : Closed Loop Simulations: Classical Design 
Approach 



The Transfer Matrix 4.6 

Transfer Matrix Compensation 14.4 

Case Study - Kidd Paper notes 

Lab - Project 2 : Open Loop Simulation, Linear MIM0 Plants 

Case Study - Kidd Paper notes 

• Multivariable Stability Introduction: 15.1, 15.2 

Liapunov’s Second Tbm. 15.3 



Lab - Project 2: Transfer Matrix Design 
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in 



11 



12 



8/19 



8/26 



9/2 



9/9 



9/16 



9/23 



Linear Systems Analysis 
Nonlinear Systems Analysis 
* Optimal Control Introduction: 



15.4 

15.6, 15.7 
7.3, 7.4, 16.1 



Lab - Project 3: Linearized Stability Analyses 



Controllability 

Observability 

Time Optimal Control 



16.2 

16.3 

16.4 



Lab - Project 3: Controllability and Observability 



Holiday 
LOR Control 

• Computer - Based Controllers: 

Introduction 

Lab - Project 3 : LOR Designs 

Sampling 

Observers 

* Model Reference, Adaptive Control 
Lab — Project 3: LOR Implementations 

• Design Summary 

♦ SOF’s 



16.5 

notes, 14.6 



notes 

notes 

16.6, 16.7 



notes 



Lab - Project 3 : Closed Loop Simulations: LOR Approach 
* No meetings, No Final 
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IE 4802 



MARINE PROPULSION CONTROL 

Project 1 - First Meek 
Classical Design: 
Multiport Analysis 
and Linearization 



12 July 1985 
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Introduction 

A typical marine propulsion emulation system is shown in the 
figure above. In this system, a water brake dynamometer is used to load 
a <*as turbine engine. The amount ot horsepower which is dissipated by t lie 
dvna^ometer is pronortion^l to the water contained in the unit and the 
speed at which the shaft rotates. The objectives of this three-week project 
are as follows: 

1. To understand how the emulator is used to mimic marine 
propulsion dynamics . 

2. To us 0 classical controller design methods to design a 
controller for the dynamometer. 

3. To evaluate the performance of the classical controller. 



2 1 



INLFT 



GAS 

GENERATOR 




N = Gas Generator Sneed 

pe; 

0 ^ = Power Turbine Toraue 
Pt 

N = Shaft Speed 



Dt 




EXHAUST 



DYNAMCHETl! 



1. In the figure above, the major components of the emulator have been 
identified for a multiport analysis. In the work for proiect 1 we will 
be considering the operation of the system at a constant water volume in 
the dynamometer. As shown in the figure, there is no gearing in this 
svstem. This means that the shaft speed is the same as the power turbine 
nT ^eed, and both of are equal to the dynamometer sneed. In this part 

of the project /you should assume that all important Inertias can be. lumped 
into the shaft component, and complete the multiport analvsis. If shaft 
speed is selected as the plant output, is the plant SISO? 



2. Construct a multiport analvsis of a machine of your choice. Comment 
on the interactions between components and overall riant function. 
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3. With the components and their inputs and outputs identified, the 
next step is to begin the modeling process. The assumption of all 
inertias being lumped into the "shaft" component allows us to use 
steady-state data for our component dynamic models elsewhere. For 
example, the steady-state performance curves for the NPS engine and 
dynamometer are attached. This data was gathered through a series of 
experiments which measured the variables in the steady-state. In our 
modeling, we will assume that this data also applies to dynamic, 
transient conditions. In this way, we can simplify the modeling task 
considerably. In this part of the project you must identify a linear 
dynamometer model of the form 

°d * f i<v- 

Recall that your design work xv r ill be for the condition of constant 
water weight (W w ) in the dynamometer. v our linearization should be 
about the operating point 

N s - 

i 

W 

v T 

4. The shaft model will be of the form 

= J + B fig , 

where J is a parameter which represents the inertia of the shaft plus 
the reflected inertias of the other components in the svstem, and 3 is 
a parameter which represents the system frictional effects. In this task 
you should identify a linearized shaft model from the above equation in 
the form 

K - <2 < "d. v- "s >■ 
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The variable 0^ is the dynamometer torque in the above equation. The 
constants in your linearization should be evaluated about the point 




N 

s 



5. The turbine torque/sneed curve which is attached is for both the gas 
generator and the power turbine, e.g. it describes the engine performance. 
In this case, a gas generator governor has been built in which regulates 
gas eenerator speed by controlling fuel flowrate (this is why N is the 

i 

input to the pas generator rather that the fuel flowrate as discussed 
in class). We will assume that the governor internal to the gas generator 
offers no significant dynamics to the plant. This assumption will allow 
you to use the turbine curves to derive a linear model of the engine of 
the form ^ 



'Pt 



f 3 ( V N » >• 



You should linearize the curves about the operating point 




N 

s 



6. Using the linear models of the plant components that vou have derived, 
find the plant transfer function which relates the plant input (gas generator 
speed) to the plant output (shaft speed). 



DYNO TORQUE (rr-LB) 

100 200 300 400 600 



DYNO TOPiQUI 



•/SPEED CUP.VE 



/ , 

V 
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ME 4802 



MARINE PROPULSION CONTROL 



I \ 



Project 1 - Second Week 



Open Loop 

Plant Model Validation 



19 July 1985 



due 



26 July 1985 
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1. Work problem 3-9-13 in Modern Control Engineering, pg 473. 



2. At the conclusion of last weeks assignment you obtained a transfer 

v 

function of the form 

n 3 a 

n bs + c 

gg 

where a, b, and c are constants which depend on J and B (the fundamental 
constants from the shaft/ inert ia model) as well as the plant operating 
point. Use the attached plant perturbation response data to estimate 
the plant inertia J and the plant friction coefficient B. Be careful to 
watch your units, n and N are in rpm while u / is in rad/sec! To check 
your work you can use the weight of the shaft (64.4 lb f , mass = 2slugs) 
along with the estimated J value in the relationship 

J = 4 m r^ . 

Thus, vour experimental value for J will allow you to calculate a 
corressponding value for the shaft radius, r. Is your computed radius 
reasonable? 

3. Previous work at the NPS has shown that the emulator operating 
curves may be roughly curve fit by the following equations 

0 d = -20 + ( (0 . 00046* (U w / lb . 6) ** 1 . 3) + 4 . OE-6) *(N g **2) 

0 . = (-725.76 + (0. 0363642*N ) ) + (0.0527 - (4 . 455E-6*N„) ) *M 

pr gg gg s 
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In this problem you should substitute the values for W^, J , and B 
into the nonlinear equations in order to find the global, nonlinear 
plant model equation of the form 



N 



s 




4. In this problem, the nonlinear equation above will be used to 
represent the "exact" dynamics of the plant and we will compare its ? 
predictions to those from the linearized dynamics developed earlier* 
The means of conducting this validation of the linear model will be 
the simulation program CSMP. It is our goal to input the same step 
sizes into each model and compare the responses of N g as a function 
of time. 



where ^^gg a "small" step. In this problem we will associate step 

smallness with the ability of the linear model to accurately predict the 
true nonlinear behavior of the plant near the operating point. In order 
to investigate this, you will need to pick appropriate values for the 



state and dynamic; a verv large step will have larger errors. Acceptable 
accuracy may be as large as + 10% in steady-state error. In our work 
we will examine the growth of the errors as the step size is increased. 



In class we discussed the relationship 



N 



gg 



N + dN 
§Si gg 




step input. A very small N will have small errors, both steady- 
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For this problem, you should submit two comparisons: one for a small 

step (steady-state error in N g about + 1%) ; and one for a step which 

goes to the limit of the operating range in N g for your dynamometer 

water weight. Either a plus or minus sten in N is OK as long as 

%% 

the same step is input to the linear and nonlinear models for 
comparison. 
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In this lab we’ll use the perturbational open-loop transfer 
function for the marine propulsion emulator to do some classical control 
design. The transfer function comes from the experimental data given 
out last week 



n 0.01 

G(s) = = 



0.2 

5 



n 



PR 



s + 

20 



s 4 - 20 



and the controller that we’d like to design is a regulator of the form 




IT' 1 * 



As users, we will input to the regulator a desired speed n 

sdes 

(note the lower case, it’s a perturbation). The controller will produce 

the plant input n^ to move the system as we desire. The objectives of 
this assignment are threefold: 

1. To design a FID controller to meet a time domain 
performance specification. 

2. To design a compensator to meet a frequency domain 

performance specification. 

3. To validate a closed loon control design using CSMP. 
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The assignment consists of three problems: 

1. What value of fieady-state controller gain, K D , is needed for a 

’’minimum system” with G(s) so that the maximum steady -state 
error is 5% ? 

2. As you observed last week, the time constant for the plant was 
very small (0,05sec.). For safety reasons, the specifications 
call for a closed-loop time constant no less than 0,5sec. To 
meet this spec you desire to slow down the closed-loop step 
response so that 67% of the final system steady-state output value 
is reached in 0,5sec, You also know that this can be accomplished 

by adding derivative control to the minimum system in the form 

\ 

rc c (s) = K p (l + t jS ) 

In this problem you should conduct an analysis to determine the 
value of Zj to give the closed-loop system time constant of 0,5sec. 

f 

You should hand in your analysis plus a CSMF run which validates 
your closed-loop design, 

3. Design a series compensator for the minimum system to give a 
bandwidth of 200rad/sec, Are stability margins a concern for this 
system? Explain your answer from a sketch of G(s) on the 
complex plane. 
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Introduction 

In this two-part lab project we will begin to look at multivariable 
control as it applies to the emulator system. The svstem shown in the 
figure above contains the dynamometer plant and three components which 
are used as control actuators. Notice that the dynamometer is separated 
into two components: the "DYNO FRICTION*’ component represents the 

dissipation of mechanical power through fluid viscous effects, and the 
"DYNO VOLUME" component represents the internal water volume of the 
unit. Notice that one-half of the dynamometer is modeled in the plant 








while the other half is modeled in the controller. Recognition of this 
effect creates a very desirable control design situation since the 
controller can thus exert change on the plant without being effected in 
return. Such high impedance connections offer an excellent place to 



separate the controller and plant for design purposes. The water flowrate 

to the dynamometer is adjusted by a controllable globe valve which is 

connected to a regulated water supply. Thus, the system inputs are 

h T , the gas generator speed, and V, the valve opening. The following 
S 

new variables appear in the figure: 






Ps 

V 



mass flowrate, Ibf/sec 
of water 

dyno water pressure, psie 
water supplv pressure, psig 
valve opening, in. 



The objectives of this two week lab are the following: 

1. To formulate multivariable system models for the components 
shown in the figure above. State variable and transfer matrix 
models will be considered. 

2. To investigate the response of the system to changes in inputs 

3. To design a controller for the system based upon the transfer 
matrix model. 

In this weeks assignment we will use a mathematical approach to 
formulating a state space model for the system, and then investigate the 
predicted response to inputs via CSMP. 
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In extending our earlier work, it can be shown that the following 

******* 

linearized model describes the dynamics of the plant 

n = 0.2 n rtrr - 20. n - 4000.0 w 

s gg s w 

Note the lower case of the variables, these are perturbations. A model 
for the dyno volume comes from the mass continuity equation 

w „ - dt - 

The NFS gas turbine dynamometer has been equipped with a regulated air 
pressure source to simplify the analysis and control design problems. 

When the air is turned on, the pressure in the dyno is = 4psig. 

The valve component can be modeled with the standard orifice flow 
equation \ 

M = K V /(P -P.) 

w v f s d' 

with P = 40psig = constant, it's also regulated 

s 

and = 1/6 lb^/sec 

This assignment Is as follows: 

1. Obtain the state equation in the standard form 
x = A x 4* B u . 

That is, what are the entries in A? and B? In order to do 

«*»* 

this, you should work in perturbational variables, e.£., v 
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(lower case) 
in C and 



is .one input while n is the other. 

| ££ 

D for the output equation 



Also find the entries 



y 


= C x + 


D 






#*► 












n 

s 


y 




w 


/ 2 a 




W w 



2. Use the INTORL statement of CSMP to integrate x = Ax + Bu 

»"* -V A. 

and plot the following: 



V 

l 


= V 1 


+ 


Y li 


= n + 

s 


N si 


Y 2 


= y 2 


+ 


Y 2i 


= W + 

w 


w , 

wi 



for = 1500 rpm 

W . = 3.5 lb. 

wi f 



for the following inputs 
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3. If N is held constant and a step input is given to V, 
ge 

under what mathematical conditions would it be possible to 
observe an oscillating response in N g ? Can a corressponding 
physical condition occur? 



I 
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This week we will conclude the plant multivariable modeline 

I 

task and design a compensator for the transfer matrix model. There are 
three parts to this assignment: 



1. h’ork problems 3-14-3 and B-14-9 in Modern Control Fneineerinp , 
pgs 713,715. 

2. Last week you found a state-space model for the XPF marine 
propulsion emulator svsten of the form 



» m - 

X 1 




n 

s 






”-20 


-4000 " 




n 

s 


# 


— 
















X 2 




V 

w 






0 

_ 


0 




w 

w 


with 




>'i 




X 1 












. V 2 




. X 2_ 






relationships 


to 


find 


the 


transfer matrix 



.2 0 

0 1 






models the emulator system above, e.^., y(s) = C (s)u(s). 



3. Design a compensator r * c (s) s0 that the operator can control eacT 
output independently , and such that the following time constants are 
realized in each control channel: 



n lag = 


3sec. 


s 




”w laB * 


lOsec 



rhat is, design so that 

y(s) = G(s)r(s) , 
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where 


V ( S ) 




n 




i 


o 




n 

sdes 




— 


s 


s: 


3s ■+• 1 












v (s’) 

o N 




w 




0 


1 

A. 








/ 




w 






10s - ij 




** ' wdes J 



Sketch the block diagram of the controlled svsten showing 
how each input channel is processed by the controller( i*e., 3 ike 
figure 14-9, pg.691 Ogata). 
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The third, and final, lab project will extend over the last five 
weeks of the quarter. The project will consist of designing an optimal 
controller for the M"" eas turbine propulsion emulator. Five steps 
will be taken to to*- final design: 



1. Develop an understanding of stability in MI' 10 systems. 

2. Conduct preliminary plant model evaluations - study 
plant controllability and observability. 

3. Perform a Linear Ouadratic Regulator (LOR) design. 

4. Determine LOR implementation strategy. 

5. Conduct LOR design validation through closed-loop 
simulation. 



This week we will* lay the foundation for understanding the LOR design 
process by investigating the stability of MIMO systems (task 1 above). In 
the following weeks, we will address the remaining steps in the list. 



This weeks assignment is composed of four parts: 



1.) Given 




’ i 


1 


3 






P = 


i 


2 


1 








3 


1 


1 




is P 

*** 


positive definite? 


(show 


your 


work) 



Also, 




* 



if V(x) 



is V(x) positive definite? 



2.) Work Droblen 3-15-2 in Modern Control Engineering. 



3 # ) Given f x. 

- x - f(x) - ; 

/" ~ ~ — i ± 

l x l + x 2_ 

use Krasovski! ! s Theorem to investigate the stability of the system 
near the origin (x^ = 0, x 0 = 0) , 

4.) Linearize the system in (3) above about the origin and compare 
the linearized conclusion with the above nonlinear conclusion. 
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This week we prepare to beein the optimal controller cesicn 



nrncess. The design recnniaue uses 2 time-domain state-space model c r 
the standard form 



X 


= 


A x 


4 - 


3 


u 


v 


= 


C x 


4* 


D 


u. 



Note that this is a linearized, nerturbational model of the plant* 

Following the identification of the state-space plant model, 
it is important to evaluate the sufficiency of the model before the 
controller is designed. Two questions about the plant model must be 
answered: Is the model controllable? and. Is is observable? Since 

the vast majority of physical plants are controllable, it is safe to 
assume that a proper model of the plant should also be controllable. 

Also, we know that the chosen state-space model is not unique for the 
given plant. Consequently, we must determine whether or not we have 
chosen a proper model form. r or controllability, we determine if 
changes in the plant inputs properly affect the plant outputs. If thev 
do, then the plant model reflects the controllability of the real 
physical plant. 

An observable plant model is one from which we can determine 
the states (x) based upon output measurements (y) . Clearly, for control 
purposes we need to *qve an adequate number and type of sensor measurements 
(y) in order to provide feedback for a regulator/controller. 
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The Durnose of this week’s assienment is to gain familiarity 
with the two preliminary plant model evaluations of controllability 
anc observability as discussed in class. The assienment has four 
parts: f 



1. Work problem 3-16-1 in Ogata. 

2. Work problem 3-16-4 in Ogata. 

3. Work problem B-16-5 in Ogata. 

4. Assess the outDut controllability and observability of the 
model of the NPS gas turbine propulsion emulator: 



w • 




- 




- 






■» — 




• 

X 1 


- 


n s 




-20 


-4000 




X 1 




• 

- X 2 . 




• 

w 




0 


0 




V 

L 4 J 


L 



“ 




" — 


n 




u. 






1 


1 













with 




Does your result agree with the test determined in 3-16-4 above? 
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Last week we studied the preliminary* evaluations of 
controllability and observabilitv of state-snace plant models. Once 
thes° determinations are made, it is nossible to begin the optimal 
controller design process. This week we will design several Linear 
Ouadratic Regulator^ r L nv c) for use in FI^O and YTYO systems. The 
assignment lias three parts: 



1. 'Cork pronlem B— 16—11 in Ogata. 



2. For the 'CPS gas turbine propulsion emulator model used 
last week, design an LOR controller based on the following 
performance index 



j 




Ox + u 




dt, 



with the following weignting matricies 



0 



*100 



o 




o 



Use the approximate methods developed in class to estimate 
the matrix. Sketch a block diagram which shows how the 

controller processes the state measurements to produce the 
plant input signals. 



3. Repeat 2 for the following weighting matricies 



*1 


0 " 


R 




“l 


0 


0 


ino 






0 


1 
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In this weeks lab we will investigate several issues of 
the control implementation design process. In this work we will assume 
that the controller has been designed already in continuous time 
and we are concerned about the effect that descret ization will have on 
system performance. The assignment has three parts: 

1, Given an e^ht bit ADC, with a maximum input signal 

of v = lOvolts, convert the input y = 6.316volts 

max 

to its corressponding digital representation. What is 
the largest error that the ADC will produce in general? 

2. Convert the following PTD controller to its approximate 
impulse transfer function form, D(z), using a backward 
difference relationship for z; 

G(s) = 10(1 + 5s + ) . 

10s 

What is the controller recursion formula for this approach? 

1. Several weeks *pm> we found the following relationship for 
the comnens itor which we designed for the N TT> S propulsion 
emulator: 

0 .6(s 4- 20 ) 6666.7 

s s 

0 0.1 

Use Euler’s method to find the approximate D^(z), and use 

D (z) to c ;nd the controller recursion formulas for each 

c 

control channel, tij and u^. 



r- (s) 

— c 
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In general, it is not always possible to measure all the 

states for an LQR state controller. This means that the control 

designer will be faced with the problem of designing an observer that 

can estimate the state values based on the available measurement 

data. In class, this problem was considered for the NFS propulsion 

emulator where the only avialable state was N , the shaft speed. The 

observer was designed to estimate both N and W based on knowledge 

s w 

of N alone (plus the plant model) . 
s 

This assignment has two parts: first (required) , repeat the 

tJPS L&fZ 

design process for the observer for the v controller designed 

in class using the value of sampling interval T = l.Osec.; and 

f4P S LOU Z 

second (optional), simulate the performance of the closed loop^system 
for the two cases of weighting matricies and the two sample intervals 
as follows, 



case 1 


0 = 


100 


0 


R = 


1 


o' 






_ 0 


1 _ 


; ~ - 

* 


0 

• 


1 







*1 


0 




P 1 


o' 


case 2 


Q = 






, R = 










0 


100 


1 ~ 


^0 


1 






w 











0.01, and l.Osec. 



0.01, and l.Osec. 



Ter the first part of the assignment, identify the entries in the 
G, JH, and matricies. Choose the observer eignvalues to be separate 

and real. For the second part of the assignment, you may find the 
attached CSMP discussion useful. The system inputs should be 



a 

27 oo 



IOOO 



'Wpesii'^ Ok*,) 



7.45 
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APPENDIX IV 
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FLUID POWER CONTROLS LABORATORY 



PROPOSED EXPERIMENTS 



I. Hydraulic Actuator Static Performance 

Goals: Determine the static performance of a 

typical rotary actuator (pump driven 
motor) . 

Tests: (a) Motor shaft locked, zero return pressure. 

Measure internal and external leakage rates 
(by collection and timing) at various supply 
pressures . 

(b) Motor shaft free and unloaded, controlled 
return pressure. Measure starting (breakout) 
pressure and running pressure at various return 
pressures . 

(c) Motor shaft free and unloaded, zero return 
pressure. Measure forward pressure vs. shaft 
speed. 

Data Reduction: Estimate internal and external leakage 

coefficients, internal friction coefficients 
(static and running) , damping coefficient. ' 
Calculate volumetric, torque, and overall 
efficiencies as function of motor speed and 
forward pressure. 

Facilities required: 

1. Instrumented flow bench and system hy- 
draulic supply. 

2. Motor (e.g. fixed displacement piston type) , 
forward and return pressure at motor ports, flow 
rate on return side, motor shaft speed. 

II. Flow control valve static performance 

Goals: Characterize a flow control valve of typical 

configuration (e.g. four-way, three-land, 
critical center or open center spool valve) . 

Tests : 



1. Valve characterization. Determine spool 
dimensions, area gradient, underlap (if any), 
etc . ■ 
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2 . 



Measure null characteristics (flow vs 
and low-side pressure drops with ?ool 



high 

centered) . 



3. 



4. 



Conduct flow versus pressure drop rests over- 
range of fixed spool displacements from full 
open to full closed. 



Compare results with theory. Estimate valve 
flow gains, flow-pressure coefficients, pressure 
sensitivies at null and show dependency on load 
pressure drop. 



Facilities required: 

1 . instrumented flow bench and system hydraulic supply. 

2 Control valve with provisions for manual spool po- 
sitioning and access to pressure measurement m 
supply and return chambers. 

3. Flow metering. Both high and low pressure capa- 
bilities are desireable. 



III. Dynamic performance of a hydraulic power element 
valve controlled piston. 
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Goals : 



To evaluate the performance of the hydraulic 
power element, determine its performance 
character istics (leakage coefficients, hydraul- 
spring rate, damping ratio, etc.) and com- 
pare these results with theory. 



Tasks : 

1. Determine necessary parameters: actuator 

volume, mass, contained volume, displacement,...) 

2. Conduct frequency response tests with no 
applied load. Determine hydraulic natural 
frequency and damping ratio at various supply 
pressures. Examine the effects of crossport 
leakage and other design modifications. 

3. Conduct stiffness (compliance) tests with 
valve stroke fixed and sinusoidal load inputs. 

4. Compare results with theory. 

Facilities required: 

1. Instrumented flow bench and system hydraulic 
power supply (constant controlled supply pressure 
at variable flow rates) . 

2. Hydraulic control valve with valve displace- 
ment controlled and measurable (valve used in 
Experiment II if suitable) . 

3. Mechanism for sinusoidal valve stroke input 
(can be home-made) . 

4. Mechanism for sinusoidal load input. 

5. Transducers for valve and load displacement, 
load force, and suitable recording/processing 
equipment . 

IV. Position Control Servomechanism 

p, 

op«rotioool I 
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Goals : 



Examine the dynamic performance of the 
servo and evaluate the effects of various 
design options. Compare observations with 
predictions . 



Tasks : 



1. Determine bandwidth and stability 
characteristics, etc. 

2. Examine the effects of under- and over- 
damping, load variations, spring loads, non- 
linearities. Induce and correct limit cycling, 
etc . 



Facilities : 

1. Flow bench with constant pressure supply. 

2. Position control servomechanism. Valve- 
controlled ram with two-stage electrohydraulic 
servovalve, position feedback, servoamplif ler . 



3. Controllable load. 

Note* These facilities should be as accessible to 
the user as possible. E.g. variable feedback gain, 
variable load damping, variable gain in pilot stage, 

etc . 



Velocity Control Servomechanism 



Velocity reference 



tQ 



Electronic 


n 


elements 




Servo 






valve 



Tachometer 






Hydraultc 

motor 






Goals: Determine system performance characteristics 

and the effects of various design options. 



Tasks-: 



Similar to experiment IV. with the addition 
of lead-lag compensation experiments and the 
effects of gear ratio at the load. 



Facilities : 



Same as in experiment IV but utilizing a rotary 
system and controlling velocity instead of position, 



FLUID POWER CONTROLS LABORATORY 



Estimated Equipment and Materials Requirements: 



Item 



Estimated Cost 



1. Fluid power bench including: 

a. Working suiface for setup and 
conduct of experiments. Includes 
mounting accomodations and excess 
fluid drainage to waste sump. 

b. Instrumented panel showing 

operating status of fluid power 
supply: Inlet and return pres- 

sures, sump temperature, return 
temperature, return flow rate. 

c. Additional pressure gages for 
indication of tapped pressures up 
to system maximum pressure. 

d. Hydraulic power supply with 
constant pressure control and 
operable at 10 gpm and 1500 psi 
supply pressure. 

e. Reservoir with capacity of 

33 Gal. Minimum and provision for 
fluid temperature control. 

f. Electric drive for above at 
230/460v, 60Hz, 3-phase. 

g. Suitable pressure connections 
for coupling experimental pressure 
to panel-mounted gages. All con- 
nections to be of low pressure-drop 
type . 



h. 



6i 



Item 



Estimated Cost 



2. Fluid power actuators suitable 

for operation at system capacity. 

a. Piston 

b. Motor 

i. Piston type 

a. Fixed Disp. 

b. Var. Disp. 

ii. Vane type. Fixed Disp. 

iii. Gear type 

3. Flow control valve - for ExpS.III, III. 
4-way, 3-land,- critical (or open) 
center spool valve with externally 
controlled spool position and inter- 
nal chambers accessible for pressure 
taps . 

4. Electrohydraulic Position Control 
Servomechanism consisting of: 

Servovalve, Actuator ' (piston) , Load 
position feedback. Items may be 
integral or separate. 

5. Electrohydraulic Velocity Control 
Servomechanism consisting of: 

Servovalve, Actuator (motor). Load 
velocity feedback. Items may be in- 
tegral or separate. 

6. Servo-amplifier suitable for application 
to items 4. and 5. above. Amplifier 
should provide for easy change of feed- 
back and feed forward gains, and se- 
lection of position or velocity con- 
trol inputs. 

7. Transducers: 

a. Direct pressure, 0-1500 psia 

b. Differential pressure, 0-800 psia. 

c. Flow- high pressure 

d. Flow- low pressure (bench mounted 
rotameter in return line) 



Item 



Estimated Cost 



8. Storage facilities: 

a. Cabinet for storing hydraulic 
hoses, connectors, etc. 

b. Cabinet for storage of valves and 
actuators . 

c. Cabinet for storage of sensors, 
transducers, and other electronic 
gear. 

9. Table working surfaces extending, in 
segments, for approximately 18 feet. 



63 



INITIAL DISTRIBUTION LIST 



No . Copies 



1. Professor Paul J. Marto, Code 69Mx 1 

Chairman, Department of Mechanical Engineering 

Naval Postgraduate School, Monterey, CA 93943-5100 

2. Associate Professor David Smith, Code 69Xh 2 

Department of Mechanical Engineering 

Naval Postgraduate School, Monterey, CA 93943-5100 

3. Professor Robert H. Nunn, Code 69Nn 2 

Department of Mechanical Engineering 

Naval Postgraduate School, Monterey, CA 93943-5100 

4. Professor Harriett Rigas, Code 62Rr 1 

Chairman, Department of Electrical and 

Computer Engineering 

Naval Postgraduate School, Monterey, CA 93943-5100 

5. CAPT George LaChance, SEA05AD 1 

Naval Sea Systems Command 

Washington, DC 20362 

6. CAPT Roger Nutting, Code 55 1 

Naval Sea Systems Command 

Washington, DC 20362 

7. CAPT Walter Ericson, Code 34 1 

Naval Engineering Curricular Officer 

Naval Postgraduate School, Monterey, CA 93943-5100 

8. LCDR Bart Everett, Code 90G 1 

Assistant for Robotics 

Naval Sea Systems Command 
Washington, DC 20362 

9. John N. Dyer, Code 06 1 

Dean of Science and Engineering 

Naval Postgraduate School, Monterey, CA 93943-5100 

10. Defense Technical Information Center 2 

Cameron Station 

Alexandria, VA 22314 

11. Dudley Knox Library, Code 0142 2 

Naval Postgraduate School, Monterey, CA 93943-5100 

12. Research Administration 1 

Code 012 

Naval Postgraduate School 
Monterey, CA 93943 



64 



DUDLEY KNOX LIBRARY - RESEARCH REPORTS 




5 6853 01057640 



U221755 



